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Abstract 

On the basis of a recently introduced generalization of the nega- 
tive binomial distribution the influence of higher-order perturbative 
QCD effects on multiplicity fluctuations are studied for deep inelastic 
e^p scattering at HERA energies. It is found that the multiplicity 
distributions measured by the HI Collaboration indicate violation of 
infinite divisibility in agreement with pQCD calculations. Attention is 
called to future experimental analysis of combinants whose nontrivial 
sign-changing oscillations are predicted using the generalized negative 
binomial law. 
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Over the past 15 years of experimental and theoretical study of multiplicity 
distributions (MDs) in a variety of collision processes a universal law emerged 
known as the negative binomial (NB) regularity [1]. The negative binomial 
distribution (NBD) proved to be very successful in describing the observed 
shape of MDs not only in full phase-space but also in restricted subdomains. 
At the same time the NBD naturally arises in branching processes which 
share many common features with the parton cascades of QCD. Thus the 
theoretical basis of the NB regularity is well founded, at least quahtatively. 

The NBD has two parameters, the shape parameter k and the mean value 
(n) . Its analytic form can be written as 

_ i r{k + n) ( e Y 

with 

Po = (1 + e)-^ and e = ^. (2) 
The probability generating function Q[z) — Z^^o-f'n-^" has the form 

G{z) = {i + {i-z)erK (3) 

A fundamental property of the distribution is that asymptotically P„ exhibits 
a nice scaling behaviour [2] 

hm P^ = l-,lj(^\ (4) 

n/{n) fixed ^ ' \ \ ' / 

which is the famous KNO scaling law [3]. The preasymptotic MDs can be 
reconstructed from the asymptotic scaling function ip{z = n/fi) via Poisson 
transform: 

Pn= r^(z)^^e-^^dz. (5) 
JO n! 

For the NBD ^(z) turns out to be a gamma distribution. In fact Eq. (5) 
is a special superposition of Poisson distributions, known as Poisson mix- 
ing in the mathematical literature [4]. The superposition is such that the 

mean values (n) of the Poisson components vary according to ip{z) which is 
called the mixing distribution. According to the theorem of Ottestad and 
Consael [4] the resulting discrete distribution P„ has probability generating 
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function equivalent to the characteristic function of ip{z) and, consequently, 
the factorial moments (cuniulants) of Pn are equivalent to the ordinary mo- 
ments(cumulants) of ■0(2;). 

During the last few years discrepancies have been found between the 
NBD and the MDs measured at the highest available energies. This is not 
surprising; it is very unlikely that a two-parameter discrete distribution can 
describe all the details of MDs in //, Ih and hh collisions. Besides the exper- 
imentally found deviations, the perturbative QCD calculations also indicate 
that the NB regularity is not general enough to account higher-order pQCD 
effects [6]. These manifest in the violation of infinite divisibility of the MDs 
in sharp contradiction to the NBD being infinitely divisible. The analysis of 
experimental data confirmed the pQCD calculations [7,8]. 

In a recent paper we have generalized the NBD by incorporating some 
pQCD-based characteristics of MDs [9]. The model involves an additional pa- 
rameter, yU, with the following physical meaning. Let us preserve the asymp- 
totic scaling form of the NBD (i.e. the gamma distribution) but not simply 
in the scaling variable z, rather, in a certain power of it, say in z'^ with /j, > 0. 
Then the asymptotic scaling function can be written as 



m 



^/^fc^M-iexp(-[A;2]^) 



(6) 



which is the generalized gamma distribution [5] . The ordinary gamma distri- 
bution is the /i = 1 special case. Obviously, the gth moment of the modified 
ip{z) involves the fractional rank q/ /i. Due to the theorem of Ottestad and 
Consael quoted before, the same rescaled rank q/ n appears in the qth fac- 
torial moment of P„ defined by the Poisson transform of the modified ip{z). 
Thus one can reproduce possible enhancement (/i < 1) or suppression (yU > 1) 
of multiplicity fluctuations with respect to the NBD by fitting the Poisson 
transform of the generalized gamma distribution to the observed P„. 

The Poisson transform of Eq. (6) can be expressed in terms of Fox's 
generahzed hypergeometric function [9]. Without going into the details we 
recapitulate the final result: 
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where 



e 



(n) m 



(9) 



and H(-) denotes the Fox- function. The probabihty generating function of 
Pn is given by 



and 



m 
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< /X < 1 
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with t = {\ — z)9. For = 1 the negative binomial distribution is recovered. 
The necessity of two separate expressions for P„ and Q{z) follows from the 
existence conditions of H(a:) discussed in [9]. 

The splitted parameter space in ji reflects an important difference between 
the two expressions for P„ given by Eqs. (7-8). According to the theorem of 
Bondesson [10] the generalized gamma density (6) is infinitely divisible only 
for < ^ < 1. For > 1 its characteristic function is entire analytic function 
of finite order and hence it must have complex zeroes which is not permit- 
ted for infinitely divisible entire characteristic functions [10]. Let us recall 
Maceda's theorem [4] which states that a discrete distribution P„ defined by 
Poisson mixing (5) is infinitely divisible if and only if the mixing distribution 
■0(2;) with z e (0, oo) is infinitely divisible. Accordingly, in case of Eq. (7) P„ 
satisfies the requirements of infinite divisibility, just as the NBD for // = 1, 
whereas the /i > 1 case given by Eq. (8) violates this feature. Thus our sim- 
ple generalization of the NBD is capable of reproducing a basic prediction 
of pQCD calculations. The additional parameter ii measures the degree of 
violation of infinite divisibility for > 1. 

To see how Eqs. (7-8) work in practice we carried out fits to the recent 
multiplicity data of the HI Collaboration measured in deep inelastic e+p 
scattering at HERA over a large kinematic region [11]. The MDs have been 
studied in pseudorapidity rj* domains of varying size in the current fragmen- 
tation region of the hadronic centre-of-mass frame. Comparison of the MDs 
to the NB and lognormal (LN) distributions showed that both models give 
an acceptable description of the uncorrected data only in the smallest pseu- 
dorapidity domains. For widening 7y*-intervals the quality of fits becomes 
progressively worse [11]. 
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In our fitting procedure numerical evaluation of the integral defining the 
Poisson transform (5) has been carried out using 96-point gaussian quadra- 
ture. It is turned out that keeping the shape parameter fixed at A; = 1 and 
fitting only (n) and fi produces reasonable description of the HI data for 
Pn tabulated in ref. [11]. Thus our theoretical P„ is the Poisson transform 
of the Weibull distribution [5]. The scale parameter of ip{z) is restricted 
to A = r(l//i)//i by the normalization condition z ■4'{z)dz = 1 [9]. The 
results of fits are collected in Table 1. As is seen the values of the ^-^e 
satisfactory, with the only exception of the fit corresponding to 1 < r^* < 3 
and 14^ = 80 ^ 115 GeV. In this case the NBD provides a better account of 
the data. Let us observe that fi > 1 and increases for widening ?7*-intervals 
signalling increasingly dominant violation of infinite divisibility. In our opin- 
ion this is the reason why the NB and LN distributions (both being infinitely 
divisible) produce progressively worse fits for large pseudorapidity domains. 
The quality of our fits is illustrated in Fig. 1 for 1 < < 4. 

In the remaining part of this Letter we consider the question of infinite di- 
visibility of MDs on the basis of combinants. According to Levy's theorem [4] 
a discrete distribution is infinitely divisible if and only if its probability gen- 
erating function can be written in the form 

g{z)=expiXigiz)-l)) (12) 

where A > and g{z) is another probability generating function. The discrete 
distributions having Q{z) of the form of Eq. (12) are known also as compound 
Poisson distributions. These can be regarded as convolutions of Poisson 
singlet, Poisson doublet, Poisson triplet, etc. distributions since Eq. (12) can 
be rewritten as 

oo 

g{z) = n cxp(C,(^« - 1)) (13) 

q=l 

i.e. as the product of the generating functions of the Poisson components 
having mean values Cq. The Cg are proportional to the probabilities given by 
g{z) with constant of proportionality A hence they can not be negative for 
infinitely divisible distributions [4]. 

After Kauffmann and Gyulassy [12] the quantities Cq acquired the name 
combinants in multiparticle phenomenology. They can be expressed in terms 
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of the probability ratios Vq = Pqj Pq according to 

Cq^Vq--Y.lC{Pq.i (14) 
1 i=l 

see [1,12,13]. Although the requirement Pq > can be a drawback for full 
phase-space analysis of MDs, in restricted subdomains the combinants have 
a number of advantages. First of all, as Eq. (14) shows, the knowledge of 
Cq requires only the finite number of probabilities Pn<q- This is extremely 
useful since testing the violation of infinite divisibility of P„ can be realized 
without the possible influence of truncating the high-multiplicity tail which 
can mimic pQCD effects [1,14]. Furthermore, one need not know the proba- 
bilities Pn themselves; the combinants can be determined directly from the 
unnormalized topological cross-sections bccaiisc they involve only ratios 
of Pn- Thus the statistical and systematic uncertainties of atot are flltered out 
by the combinants. They also share some common features with the factorial 
cumulants, e.g. Ci = (n) and Cq>2 — for the Poisson distribution, further, 
both quantities should be nonnegative for infinitely divisible distributions. 

Since the generalization of the NBD obtained by the Poisson transform of 
Eq. (6) can violate infinite divisibility it is of interest to study the behaviour of 
its combinants. For < 1 one expects positive Cq in all orders but for > 1 
sign-changing oscillations may appear in the g-dependence of Cq similarly 
to the factorial cumulant-to- moment ratios Hq of the distribution [9]. We 
have calculated the combinants of the generalized NBD with the help of 
Eq. (14) using the same numerical integration procedure for the evaluation 
of Pn as described earlier. A typical result for the behaviour of log \Cq\ over 
the /i-q plane is shown in Fig. 2a for k = 3/2 and (n) = 10 (with the same k 
the behaviour of Hq is studied in [9]). The peculiar structure for > 1 is due 
to sign-changing oscillations of Cq. At a fixed q the neighbouring bumps are 
Cq intervals of opposite sign in /x. For /i < 1 Cg is always positive as expected. 
In Fig. 2b slices are shown with g = 5, 20 (left) and with fi = 1, A (right). It 
is seen that the pattern of oscillations is nontrivial, i.e. not alternating in sign 
as q takes even/odd values. Qualitatively similar sign-changing oscillations 
of Cq occur for a different choice of (n) and k. 

Finally we provide a further example for the usefulness of combinants 
in the analysis of MDs. It is often stated that the factorial cumulant-to- 
moments raios Hq are very sensitive to tiny details of the high-multiplicity 
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tail of Pn (and thus to truncation effects as well). Let us consider this 
widespread opinion somewhat closer for the NBD. The unnormalized fac- 
torial moments of the NBD are given by 

where is given in Eq. (2). The unnormalized factorial cumulants fq take 
the form 

fq = kV{q)e'' (16) 
hence for Hq = fq/^q one obtains [6] 

Hq = k'^^^^kBik,q) (17) 

with B(-) denoting the Euler beta-function. The combinants of the NBD are 
given by [15,16] 

Comparing Eqs. (1) and (15) we get 

' ' r(k) \i + e) (! + «)« ^ ' 

whereas the comparison of Eqs. (18) and (16) yields 

= (20) 

Thus one arrives at the rather unexpected result 

C 1 '^"^ P • 

H. = ^ = 1-IT.^C.^ (21) 
' q ^ i=i 

which shows that the factorial cumulant-to- moment ratios Hq of the NBD are 
completely insensitive to the Pn>q tail of the distribution and carry essentially 
the same information as the combinants. This is a very special property of 
the NBD which does not hold in general for infinitely divisible distributions. 
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In summary, we have investigated the violation of infinite divisibihty of 
MDs. This particular feature plays a central role in recent studies of mul- 
tiplicity fluctuations since the higher-order pQCD calculations predict its 
appearance for MDs measured in full phase-space and in restricted subdo- 
mains. As a consequence, departures arise from the NB regularity. 

In a previous paper we have generalized the NBD by extending the valid- 
ity of its asymptotic scaling form to a certain power of the scaling variable. 
By this modification one can reproduce possible suppression/enhancement 
of multiplicity fiuctuations with respect to the NBD as well as the violation 
of infinite divisibility of MDs. Fitting the generalized NBD to the recent 
experimental data for P„ measured by the HI Collaboration in deep inelastic 
e'^p scattering at HERA we have obtained good agreement. According to 
our results the violation of infinite divisibility of P„ is increasingly dominant 
in widening pseudorapidity intervals. 

Investigating the combinants of the generalized NBD we have found non- 
trivial sign-changing oscillations of these quantities for the subdomain of 
parameter space violating infinite divisibility. In our view the higher-order 
pQCD effects on multiphcity fluctuations can be studied most effectively 
by measuring the combinants in restricted phase-space volumes. The con- 
straints of conservation laws are less pronounced, hence the dynamical ef- 
fects are more visible, and the sign-changing oscillations of the combinants 
are not corrupted by finite statistics effects. The pQCD prediction for the 
^-dependence of Cq (location of the first minimum, etc.) would be very im- 
portant to identify properly the dynamics underlying the oscillations. 
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77*-interval 


W (GeV) 


/i 


in) 


XV d.o.f. 


l<r]* <2 


80 115 


1.96 ±0.08 


2.44 ±0.04 


18/12 


115^ 150 


1.82 ±0.07 


2.50 ±0.04 


5/13 


150 ^ 185 


1.87 ±0.08 


2.61 ± 0.05 


12/13 


185 ^ 220 


1.84 ±0.09 


2.64 ±0.06 


7/13 


1 < 77* < 3 


80^115 


2.46 ±0.07 


4.88 ±0.06 


44/16 


115^ 150 


2.09 ±0.06 


5.05 ±0.07 


25/17 


150 ^ 185 


2.07 ±0.07 


5.29 ±0.08 


9/19 


185 ^ 220 


2.10 ±0.08 


5.33 ±0.09 


7/20 


1< 77* < 4 


80 ^ 115 


3.70 ±0.11 


6.12 ± 0.06 


17/17 


115 150 


3.18±0.11 


7.02 ±0.07 


18/19 


150 ^ 185 


2.92 ±0.09 


7.49 ± 0.08 


16/21 


185 ^ 220 


2.75 ±0.10 


7.67 ±0.09 


12/21 


1 < 7^* < 5 


80 115 


4.66 ±0.23 


6.87 ±0.06 


12/17 


115^ 150 


4.44 ±0.18 


7.70 ± 0.06 


26/20 


150 ^ 185 


4.16 ±0.16 


8.39 ±0.08 


19/21 


185 ^ 220 


4.08 ±0.18 


8.79 ± 0.08 


16/22 



Table 1. Results of fits to the HI data for P„ with the Poisson transform 
of the generalized gamma distribution (6) for fixed shape parameter k — 1 
(WeibuU case). The errors are only statistical. 



FIGURE CAPTIONS 

Fig. 1: The best-fit theoretical distributions for 1 < 77* < 4, see the text 
for details. The corresponding parameters are collected in Table 1. The 
displayed errors are only statistical. 

Fig. 2a: Sign-changing oscillations of the combinants Cg of the Poisson 
transformed generalized gamma distribution. The parameters kept fixed 
are k — 3/2 and (71) = 10. For clarity only the odd-rank combinants are 
displayed. 

Fig. 2b: Slices through Fig. 2a with g = 5, 20 (left) and with /i = 1, 4 (right). 
The fjL — l case (smooth curve) corresponds to the NBD. 




Figure 2b 



